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Abstract 

We present measurements of the branching fractions of three-prong and five-prong r decay modes 
using a sample of 430 million r lepton pairs, corresponding to an integrated luminosity of 468 fb _1 , 
collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ storage rings at SLAC. 
The t~ — > (37r)~ r\v T , t~ — > (37r)~ uju t , and t~ — > n~ fi (1285)f T branching fractions are presented, 
as well as a new limit on the branching fraction of the second-class current decay t~ — > ir~ rj (958)^ T . 
We search for the decay mode t~ — > K~ r\ (958)f T and for five-prong decay modes with kaons, and 
place the first upper limits on their branching fractions. 

PACS numbers: 13.35.Dx, 14.60.Fg 



I. INTRODUCTION 

Study of the three-prong and five-prong decay modes 
of the r lepton, where "prong" refers to the number of 
charged hadrons (it or K) in the final state, allows one to 
test the Standard Model and search for evidence of new 
physics. The large r lepton data sample collected by the 
BABAR experiment provides an opportunity to perform 
a comprehensive study of rare, high multiplicity decay 
modes and to search for forbidden processes. 

We present measurements of the t~ — > (3n)~ t]v T , 
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t~ — > (3ir)~uiis T , and t~ — > ir~ f\V T branching fractions. 
We use the primary decay modes of the 77, a; (782), and 
/i(1258) mesons: r/ — > 77, 77 — > 7T + 7r~7r°, 77 — > 3ir°, 
uj — > 7T~7r + 7r , fi — > 27T + 27r~, and fx — > Tr + ir~i]. No 
other narrow resonances are observed. We measure the 
branching fractions of the non-resonant decays, where 
the non-resonant category includes possible contributions 
from broad resonances. We present a new limit on 
the branching fractions of the second-class current de- 
cay r~ — > 7r~?/(958)^ r , and the first limits on the al- 
lowed first-class current decays r~ — > K~rj' (958) v T and 
r~ — > 7r~7r°?7'(958)^ T . Finally, we present the first limits 
on the branching fractions of five-prong decay modes in 
which one or more of the charged hadrons is a charged 
kaon. Note that the branching fractions exclude the con- 
tribution of K° s — > tt + tt~ decays. Throughout this paper, 
charge conjugation is implied. 

This analysis is based on data recorded with the BABAR 
detector at the PEP-II asymmetric-energy e + e~ stor- 
age rings operated at the SLAC National Accelerator 
Laboratory. With an integrated luminosity (£) of 424 
and 44 fb^ 1 recorded at center-of-mass (CM) energies of 
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10.58 GeV and 10.54 GeV, respectively, and an e + e~ — >• 
t + t~ cross section of a T + T - = (0.919 ± 0.003) nb [l|, the 
data sample contains 430 million r pairs. 

The BaBAR detector is described in detail in Ref. p|. 
Charged-particle momenta are measured with a five-layer 
double-sided silicon vertex tracker and a 40-layer drift 
chamber, both operating in the 1.5 T magnetic field of a 
superconducting solenoid. Information from a detector of 
internally reflected Cerenkov light is used in conjunction 
with specific energy loss measurements from the tracking 
detectors to identify charged pions and kaons Q. Pho- 
tons are reconstructed from energy clusters deposited in a 
CsI(TI) electromagnetic calorimeter. Electrons are iden- 
tified by combining tracking and calorimeter information. 
An instrumented magnetic flux return is used to identify 
muons. 

The background contamination and selection efficien- 
cies are determined using Monte Carlo simulation. The 
r-pair production is simulated with the KK2F event gen- 
erator |4J. The t decays, continuum qq events (q = udsc), 
and final-state radiative effects are modeled with the 
Tauola @, JETSET @, and Photos Q generators, re- 
spectively. Dedicated samples of t + t~ events are cre- 
ated using the Tauola or EVTGEN || programs, with 
one of the r leptons allowed to decay to any mode while 
the other r decays to a specific final state. The detec- 
tor response is simulated with GEANT4 j|. All Monte 
Carlo events are processed through a full simulation of 
the BABAR detector and are reconstructed in the same 
manner as the data. 



II. EVENT SELECTION 

The t pair is produced back-to-back in the e + e~ CM 
frame. The decay products of the two r leptons can 
thus be separated from each other by dividing the event 
into two hemispheres: the "signal" hemisphere and the 
"tag" hemisphere. The separation is performed using 
the event thrust axis which is calculated using all 
charged particle and photon candidates in the event. 

We select events where one hemisphere (tag) contains 
exactly one track while the other hemisphere (signal) con- 
tains exactly three or five tracks with total charge oppo- 
site to that of the tag hemisphere. The event is rejected 
if any pair of oppositely charged tracks is consistent with 
a photon conversion. The component of the momentum 
transverse to the beam axis for each of the tracks must 
be greater than O.lGeV/c in the laboratory frame. All 
tracks are required to have a point of closest approach 
to the interaction region less than 1.5 cm in the plane 
transverse to the beam axis and less than 2.5 cm in the 
direction along that axis. This requirement eliminates 
mesons that decay to 7r + 7r~ at points distant from 
the e + e~ collision point. 

To reduce backgrounds from non- r-pair events, we re- 
quire that the momentum of the charged particle in the 



tag hemisphere be less than 4 GeV/c in the CM frame and 
that the charged particle be identified as an electron or 
a muon. The qq background is suppressed by requiring 
that there be at most one energetic (E > 1 GeV) electro- 
magnetic calorimeter cluster in the tag hemisphere that 
is not associated with a track. Additional background 
suppression is achieved by requiring the magnitude of 
the event thrust to lie between 0.92 and 0.99. 

Neutral pion and eta candidates are reconstructed from 
two photon candidates, each with energy greater than 
30 MeV in the laboratory frame; the invariant mass of 
the 7T° (77) is required to be between 0.115 (0.35) and 
0.150 (0.70) GeV/c 2 . Neutral pion candidates are recon- 
structed in the signal hemisphere. If a photon candi- 
date meets the invariant mass requirement with multi- 
ple photon candidates, then the neutral pion candidate 
with invariant mass closest to the nominal ir° mass [TTj 
is selected. The search for additional pion candidates 
is repeated using the remaining photon candidates. The 
residual photon clusters in the signal hemisphere are used 
to search for 77 —> 77 candidates. In the case of multiple 
t] candidates, the candidate with invariant mass closest 
to the nominal 77 mass is selected. We reject events in 
which the invariant mass M formed from the system of 
charged particles, 7r°, and 77 candidates, all in the signal 
hemisphere, exceeds 1.8GeV/c 2 . 

The branching fractions are calculated using the ex- 
pression B = Nx/(2Ne) where Nx is the number of can- 
didates after background subtraction, N is the number 
of r pairs produced, and e is the selection efficiency. N is 
determined from the product of the integrated luminos- 
ity and the e + e~ — > t + t~ cross section. The uncertainty 
of N is estimated to be 1%. The selection efficiencies are 
determined from the signal Monte Carlo samples. The 
uncertainty on the selection efficiencies includes 0.5% per 
track on the track reconstruction efficiency, as well as par- 
ticle identification (PID) selection uncertainties. From 
studies conducted on real and simulated events, the un- 
certainties on the charged particle identification selectors 
are estimated to be 1% for electrons, 2.5% for muons, 
0.5% for pions, and 1.8% for kaons. The combined elec- 
tron and muon particle identification uncertainty is esti- 
mated to be 1.6% based on the composition of the event 
samples. The uncertainty on the 7r° — > 77 and 77 — > 77 
reconstruction efficiency is estimated to be 3% per can- 
didate. 



III. RESULTS 

We present measurements of r decays to a system with 
77, fx and lu resonances in Sections A, B, and C, respec- 
tively. Decays with these resonances do not account for 
all three-prong or five-prong r decay modes, as discussed 
below, and we present measurements of the r branching 
fractions through non-resonant modes in Section D. Sec- 
tion E presents a search for r decays with an 77' (958) 
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FIG. 1: The 77, 7r + 7r _ 7r°, and 3tv° invariant mass distribu- 
tions for r~ — > 2ty~ -k + -qv T decay candidates, and the 7r + 7r _ 7r° 
invariant mass distribution for t~ — > -K~2-K°riv T decay candi- 
dates, after all selection criteria are applied. The solid lines 
represent the fit to the r\ peak and background. The dashed 
lines show the extrapolation of the background function under 
the 7] peak. 



meson, while Section F presents a search for decays with 
cither one or two charged kaons. 



A. t — > (371") r\v T 

The t~ — ¥ 2Tr~n + r)i' T mode is studied in the rj— > 77, 
Tf—¥ 7r + 7T~7r , and 77 — >• 3ir° final states, while the r~ — > 
7r _ 27r°?7^ T mode is studied in the 77— > 7r + 7r~7r final state. 

The event yields arc determined by fitting the 77 mass 
peak in the 77, 7r + 7r~7r°, and 3ir° invariant mass distri- 
butions (see Fig. []}. The fit uses a Novosibirsk function 
12] (a Gaussian distribution with a tail parameter) for 
the 77 and a polynomial function for the background. 

The Monte Carlo simulation indicates that some of 
the entries in the rj peak are from e + e~ — > qq events. 



Control samples, obtained by reversing the requirement 
on the invariant mass of the observed decay products 
(M > 1.8GeV/c 2 ), are used to validate the background 
estimate. The expected background is corrected by the 
ratio of data to Monte Carlo events, and the statistical 
uncertainty of the ratio is included in the background 
systematic uncertainty. This method of validating the 
qq background estimate is used for all decays and is not 
mentioned in the later sections. 

The reconstruction efficiencies are determined from fits 
to the signal Monte Carlo samples. The t" — > -K~2ifir\v T 
sample is generated using a phase-space model for the 
final-state particles. The r~ — > 2ix~ "7r + r\v r sample is 
composed of t~ — > ir~ f\v r (fi — > 7r + 7r — 77) decays and 
decays without an intermediate resonance. The r~ — > 



2n Tr + rjv T (excluding /1) and r 



f\v T efficiencies 



are the same for 77 — > 7r + 7r~7r u and 77 — jv 3ir° events, 
whereas a slight difference is observed for 77 — s> 77 events 
and a 2.5% uncertainty is added to the selection efficiency 
systematic for this mode. In addition, a 4% uncertainty 
is added to the t~ — > 2n~Tr + riv T selection efficiency for 
the 77 — > 7T + 7r~7r° mode to take into account variations 
observed for different fits. 

The three determinations of the t~ — > 2~k~ ~tt^ r\v T 
branching fraction are found to be in good agreement (see 
Table HI and we therefore calculate a weighted average. 
The statistical and systematic uncertainties on the aver- 
age are obtained by combining the individual uncertain- 
ties in quadrature, accounting for correlations between 
the systematic terms. The weighted average (inclusive of 
t~ — > n~ f±v T ) is found to be 

B{r~ -> 2ii-Tt + r,v T ) = (2.25 ± 0.07 ± 0.12) x 10~ 4 . 

Hereinafter, when two uncertainties are quoted, the first 
is statistical and the second is systematic. The average 
branching fraction (exclusive of t~ — > 7r~ fiv T ) is deter- 
mined to be (0.99 ± 0.09 ± 0.13) x 10~ 4 and is obtained 
using the branching fraction (inclusive of t~ — > n~ f\v T ), 
given above, and subtracting the product branching frac- 
tion B(t~ — ¥ 7r~ f\v T ) x B(fi — > TT + TT~ri) presented in the 
next section. 

The r~ — > n~2iJpr\v T branching fraction is found to be 



"2ttVt) = (2.01 ±0.34 ±0.22) x 10" 



Naively, the ratio of the — > 2n~ ii + r\v T to r~ — > 
'K~2'K°rjv T branching fractions is expected to be two if 
the decay is dominated by the r~ — > tt~JiV t decay mode 
(based on the /1 branching fractions ill,]). The data do 
not support this expectation. 

Our previous measurement of the t~ — » 2TT~TT +r qv T 
branching fraction (1.60 ± 0.05 ± 0.11) x 10~ 4 [II], which 
is based on the 77— > 77 mode only, is superseded by this 
measurement. The fit used in the previous analysis was 
performed using a narrower range in the invariant mass 
distribution, and a background parametrization that is 
less accurate than in the present study. 

The r~ — s> 2Tx^-K + rjv T and r~ — > n~2iJpr\v T branching 
fractions are in good agreement with the results from 
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TABLE I: Results and branching fractions for r — > (3n) r\v T decays. 



Decay Mode 


T — )■ 27T 7V~^T}h / T 


T — r 2"7T 7V~^T)h / T 


T — ^ 27T 7T~^TjL' T 


T — r TV 27V^7]U T 






TJ — y TV 7T 7T ^ 


Tj — S> 37T° 


TJ — y 7V~^~ 7Y 7T^ 






2S7 + 12 + 18 


2 54 + 27 + 25 


2 01 + S4 + 22 


Data events 


2887 ± 103 


1440 ± 68 


315 ± 34 


381 ± 45 




107/76 


60 /52 


31/34 


95/75 


Selection efficiency 


(3.83 ± 0.11)% 


(2.97 ± 0.12)% 


(0.42 ± 0.01)% 


(0.75 ± 0.02)% 


RapVorni inn p-iT-p-ri-f q 


131 ± 29 


65 ± 38 


13 ± 7 


83 ± 12 


Svcj't'PTnafip mi pprt ft infipQ [%i 










TVflplcin cr pffipiptipv 

1 1, QiL' Iv l XL Li GlllLlCllv y 


2.7 


3.8 


2.7 


2.7 


7T° and r\ PID 


3.0 


3.0 


9.0 


9.0 


Pion PID 


1.5 


2.5 


1.5 


1.5 


Lepton-tag PID 


1.6 


1.6 


1.6 


1.6 


TV 


1.0 


1.0 


1.0 


1.0 


Selection efficiency 


3.9 


4.0 


2.8 


2.7 


Background 


1.0 


2.8 


2.3 


4.0 


B(J7-» 77) 


1.0 








B(?7->- Tr+TT^TT ) 




1.8 




1.8 


/3(r?->- 3tt°) 






0.9 




Total (%) 


6.3 


7.4 


10 


11 



the CLEO Collaboration, (2.3 ± 0.5) x 1CT 4 and (1.5 ± 
0.5) x 10~ 4 , respectively [Tj]. Li predicts a larger t~ — > 
2-K~7r + rjv T branching fraction, 2.93 x 10 -4 [i~5| . 



B. 



TV f\V T 



The branching fraction of — > n~ j\v r and the mass 
of the /i meson are measured using the fx — > 2~K + 2 r K~ and 
fx — > 7r + 7r~?7 decay modes, where the fx — > ir + T:^r] decay 
is reconstructed using 77— s> 77, 77 — > 7r + 7r~7r°, and 77^ 37r° 
events. The criteria used to select the r~ — > tt~ fxv T 
decays for the branching fraction measurement are de- 
scribed earlier. We modify the selection for the mass 
measurement, dropping the requirement that the track 
in the tag hemisphere be a lepton and the restriction on 
the number of photon candidates in the tag hemisphere, 
to increase the size of the event sample. 

The numbers of r~ — > n~ f\V T candidates are deter- 
mined by fitting the fx peak in the 27r + 27r~ and 7r + 7r~?7 
invariant mass distributions (see Fig. [2]). The /1 line- 
shape is expected to be a Breit-Wigner distribution, mod- 
ified by the limited phase space. Previous studies show 
that the /1 —> a^~7r + (a^ (980) — > n~r]) channel appears 
to account for all fx — > tt + tt~i] decays [HI]. The mass of 
the 7rao (980) system and the r mass provide a lower and 
upper limit, respectively, on the /1 lineshape. We use 
the four- vectors of the charged pion and do (980) from 
the EVTGEN generator to determine the simulated /1 
lineshape and find it to be a close approximation to the 
Breit-Wigner expectation. The fx peak is fit using this 



lineshape convolved with a Gaussian distribution to take 
into account the effects of the detector resolution. The 
results of the fits are presented in Table HU There is no 
evidence for peaking background from qq events or other 
t decays. 

The product of the r~ — > tt~ f\v T and /1 — > 2tt + 2it~ 
branching fractions, and the product of the t" — > n~ f\v T 
and /1 — > 7r + Tr~r] branching fractions, are measured to be 

B(r~-> tt-/ii/ t )B(/i-> 2^+2tt-) 
= (5.20 ±0.31 ±0.37) x 10~ 5 , 
B{t" -» -K~fiV T )B{fi ->■ Tr+Tr^ry) 
= (1.26 ±0.06 ±0.06) x 10~ 4 , 

respectively, where the second result is the weighted aver- 
age of the three rj modes. The B(t~ —> n~ f\v T ) branch- 
ing fraction is determined to be (4.73±0.28±0.45) x 10~ 4 
and (3.60 ± 0.18 ± 0.23) x 10~ 4 , as obtained by dividing 
the product branching fractions by £>(/i — > 27r + 27r~) = 
Q.110+°;™6 and -> = 0.349t°;°^ ^ re - 

spectively. 

Our two measured values for the r _ — > n~ j\v T branch- 
ing fraction are consistent with each other to within two 
standard deviations of the combined statistical and sys- 
tematic uncertainties. The ratio of the product branch- 
ing fractions is used to determine the ratio of the fx — > 
2tt + 2tt~ and fx — > 7r + 7r~?7 branching fractions as 



27T+27T-) 
-)• 7T7T7/) 



0.28 ±0.02 ±0.02, 
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TABLE II: Results and branching fractions for r — > n f\v T decays. 



Decay Mode 


27T+271-- 


/l — > 7i~ + 7T 77 


/1 — > n + TT 77 


/l — > TV + TT 77 






77— >• 77 


77— > 7!" + 7r _ 7T 




Branching fractions (10 -4 ) 










B(t~-» Tr~fiv T )B(fi -> 2-k + 2%-) 


0.520 ±0.031 ±0.037 








B(t~ — > 7r~/lI/ T )B(/l — > 7r _ 7T + 77) 




1.25 ±0.08 ±0.07 


1.26 ±0.11 ±0.08 


1.33 ±0.39 ±0.20 


Data events 


3722 ± 222 


1605 ± 94 


731 ± 62 


197 ± 59 


X 2 /NDF 


77/62 


50/43 


61/55 


39/43 


Selection efficiency 


(8.3 ±0.1)% 


(3.75 ± 0.04)% 


(2.97 ±0.05)% 


(0.53 ± 0.06)% 


Systematic uncertainties (%) 










Tracking efficiency 


3.8 


2.7 


3.8 


2.7 


7T° and 77 PID 


- 


3.0 


3.0 


9.0 


Pion PID 


2.5 


1.5 


2.5 


1.5 


Lepton-tag PID 


1.6 


1.6 


1.6 


1.6 


N 


1.0 


1.0 


1.0 


1.0 


Selection efficiency 


0.6 


1.1 


1.6 


11 


Fit model 


5.0 


2.7 






/3(t? ^ 77) 




0.7 






B{r)^ n+TY-n ) 






1.2 




B(t]^ 3tt°) 








0.9 


Total (%) 


7.0 


5.6 


6.1 


15 



TABLE III: Results of fits for the mass of the /1 resonance in r — > 7r fxv T decays. The errors are statistical. 



Decay Mode 


Monte Carlo 


Data 


Data 




(generator - fit) 


(fit) 


(corrected) 




(GeV/c 2 ) 


(GeV/c 2 ) 


(GeV/c 2 ) 


27T+27T" 


0.00074 ± 0.00008 


1.28031 ± 0.00067 


1.28105 ± 0.00067 


fl — > -K + 1T~ T] 








77— > 77 


0.00292 ± 0.00040 


1.27775 ± 0.00045 


1.28067 ± 0.00060 


77— > 7r + 7T _ 7T 


0.00018 ± 0.00020 


1.27787 ± 0.00080 


1.27805 ± 0.00082 


7/^ 37T° 


0.00347 ± 0.00033 


1.28036 ± 0.00335 


1.28383 ± 0.00337 



where B(fi — > ttttt]) = 1.5 X £>(/i —> Tt + Tt~ri) based 
on isospin symmetry. This agrees with average value of 
0.41 ± 0.14 quoted by the Particle Data Gr oup [Tl| but 
disagrees with their fit value of 0.63 ± 0.06 11] . 

The systematic uncertainties of the branching fractions 
are listed in Table [TTJ We observe that the number of 
events in the /1 peak in the fx — > 2tt + 2tt~ sample varies 
by 5% for different background shapes. This variation is 
included as a systematic uncertainty. We also observe 
that the selection efficiency obtained from the Monte 
Carlo simulation exhibits a slight dependence on whether 
the f\ decays via the /1 — > a^n + or the /1 — > tt + tt~i] 
mode, and the variation is included as a systematic un- 
certainty (listed under "Fit model" in Table HI)) . 

The t~ — > tt~ f\V T branching fraction using the /1 — > 




2tt + 2tt~ mode is consistent with the previous BABAR 
measurement, which is also based on the f\ — > 2tt + 2tt~ 
mode [18(. CLEO published a branching fraction of 
(5.8^3 ± 1.8) x 10~ 4 |19j and Li predicts a branching 
fraction of 2.9 x 10~ 

The fi mass is determined by fitting the peak with a 
non-relativistic Breit-Wigner function, which was used in 
previous measurements of the /1 mass [llj. As a cross 
check, we use the energy-momentum four-vectors from 
the generator Monte Carlo simulation, and we find the 
fitted mass value to be consistent with the input mass 
value. 

We fit the invariant mass distribution in the fully re- 
constructed Monte Carlo samples to determine whether 
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FIG. 2: The 2ir 2ir~ and tt + it~t] invariant mass distribu- 
tions for t~ — > 2ir~ -k + r]v T decay candidates after all selection 
criteria are applied. The lower three plots are for the r\— > 77, 
77— > ty + tt~ 7T° , and r\— > 37r° decays. The solid lines represent 
the fit to the /i(1285) peak and background. The dashed lines 
show the extrapolation of the background function under the 
/1 peak. 



the result of the fit differs from the input mass of the 
Monte Carlo generator. The difference is used to correct 
the value of the invariant mass of each channel obtained 
from the fit and the uncertainty in the difference is in- 
cluded as a systematic error. 

Table Mil and Fig. [3] show the results of the fits to 
the data. The last column of the table presents the 
mass after the application of the reconstruction cor- 
rection factor. The average of these results is Mf 1 — 
(1.28025 ±0.00039) GeV/c 2 , where the error is statistical. 

Previous BABAR analyses have measured the invariant 
mass of resonances to be approximately 1 MeV/c 2 less 
than the values of Ref. [ll|. This shift is observed in 
the measurement of the mass of the /1 meson [2l| and 
the r lepton (22j. The shift is attributed to the abso- 
lute energy and momentum calibration of the detector. 



> 71 71 T] 



T| — > 7l + 7["7[ u 

n — > y y 




1.26 1.265 1.27 1.275 1.28 1.285 1.29 

Mass (GeV/c 2 ) 



FIG. 3: Compilation of our measurements of the /1 mass. 
The solid line is the weighted average and the shaded area is 
the one-standard-deviation region. 
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FIG. 4: The fits to the iv peak in the -k tt~-k invariant mass 
distributions for t~ — > 2n~ -k + uv T and r~ — > ■k~2-k°ujv t de- 
cay candidates after all selection criteria are applied. The 
solid lines represent the fit to the lo peak and background. 
The dashed lines show the extrapolation of the background 
function under the ui peak. 



We measure the calibration correction factor by fitting 
the rj, lj, rj ', D° and D*~ states using data samples that 
have one track in the tag hemisphere and either three or 
five tracks in the signal hemisphere. No other selection 
criteria are applied. The peak masses are found to be 
lower than the known values by (—0.91 ± 0.10) MeV/c 2 
and the values are independent of the mass of the reso- 
nance. The calibration correction factor is applied to the 
invariant mass and its error is included in the systematic 
uncertainty. 
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We determine the mass of the /i(1258) meson to be 

M h = (1.28116 ± 0.00039 ± 0.00045) GeV/c 2 . 

The systematic uncertainty includes the reconstruction 
uncertainty and the calibration uncertainty. This result 
is in good agreement with the value (1.2818 ± 0.0006) 
GeV/c 2 inRef. [HI- 
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FIG. 5: The 3ir° , -k + -k~tt°, and 2n~ ir + 3ir° invariant mass 
distributions in r~ — > 2tt^tt + 3tv°u t decay candidates. The 
predictions of the Monte Carlo simulation are shown for the 
resonant and non-resonant r decays, and the background from 
other r decays and qq events. The resonant decays include 
decays with the correct topology and a resonance (77, /1 or oj) 
in the final state. 



C. 



(37r) U)V t 



We measure the t~ — > 2n~ir + (jjv T and r~ — > ix~2'K i ujv T 
branching fractions. The number of events is determined 
by fitting the to peak in the 7r + 7r~7r° invariant mass dis- 
tributions (see Fig. |4]) with a Breit-Wigner distribution, 
which is convolved with a Gaussian distribution to take 
into account the detector resolution. The resolution pa- 
rameter of the Gaussian distribution is determined using 
a data control sample consisting of qq events, and is fixed 
in the fit. A polynomial function is used to fit the back- 
ground. The results are presented in Table HVl 

Approximately 10% of the events in the t~ — > 
2tt~ 7t + iqv t channel arc backgrounds from other r de- 
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Non-resonant 

Resonant 
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1 1.2 1.4 
7t + tT Mass (GeV/c 2 ) 




1.2 1.4 1.6 1.8 2 
37T 27t+ Mass (GeV/c 2 ) 



FIG. 6: The tt + tt~ , 2tt + 2tt~ , and ?>n~2ii + invariant mass 
distributions in r~ — ¥ 3Tv^2-K + f T decays. The predictions 
of the Monte Carlo simulation are shown for the resonant 
and non-resonant r decays, and the background from other 
r decays and qq events. The resonant decays include decays 
with the correct topology and a resonance (77, /1 or ui) in 
the final state. The non-resonant decays are generated using 
r~ — > a~[v T events. The differences between the data and 
Monte Carlo predictions are discussed in the text. 



events. The backgrounds are subtracted before calculat- 
ing the branching fraction. 

The t~ — > 7T~27r°ajjy r sample has substantial contribu- 
tions from t~ — > it~ojv t and t~ — > n~ir°ujv T decays. The 
background is estimated with the Monte Carlo simulation 
and verified using data and simulation control samples. 
The control samples follow the nominal selection criteria 
but select one or two 7r° instead of three ir° mesons. 

The branching fractions are found to be 



(8.4 ±0.4 ±0.6) x 10~ 5 , 
(7.3 ±1.2 ±1.2) x 10- 5 . 



cays (primarily r 



7T 7T v Wiy, 



decays) and e" 1 



qq 



B(t -} 2ir it + ijjv t ) 



The systematic uncertainties are listed in Tabic lIVl 

The r~ — > 2i:~Tt + ujv T and t~ — > tt~2it uv r branch- 
ing fractions are consistent with the results from CLEO, 
(1.2 ± 0.2 ± 0.1) x 10~ 4 and (1.4 ± 0.4 ± 0.3) x 10~ 4 , 
respectively [l4j]- Gao and Li suggest that this mode is 
dominated by the (irpuj) state and predict a branching 
fraction in the range of 1.8 - 2.1 x 10" 4 with the two 
modes (t~ — > 2ir~ir + Luv T and r~ — > tt~2'k (jJv t ) having 
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TABLE IV: Results and branching fractions for r — > (3n) uiu T decays. 



Decay Mode 


T — ^ 27T 7T^~L<Ji^ r 

- + o 
— ^. 7f 7T 7T 


T — Y TV ln^UJVr 

- + o 

CO — ¥ 7V TV TV 


Rrflnptiino' frfictinvm fid J ^ 




7 S + 1 2 + 1 2 


Data events 


2372 ± 94 


1135 ± 70 


v 2 / /V D F 


55/44 


42/44 


Selection efficiency 


(3.27 ± 0.03)% 


(0.75 ± 0.01)% 


XJdCrLgl U LII1L1 


9^7 4- 71 




Systematic uncertainties (%) 






Tracking efficiency 


3.8 


2.7 


7T° and T] PfD 


3.0 


9.0 


Pion PID 


2.5 


1.5 


Lepton-tag PfD 


1.6 


1.6 


AT 


1.0 


1.0 


Selection efficiency 


0.8 


1.8 


Background 


3.4 


14 




0.8 


0.8 


Total (%) 


6.8 


17 




FIG. 7: The -k + -k~ , 7r + 7r~7r°, 27r~27r + 7r°, and 2>n~ 2-k + n invariant mass distributions in t~ — > 37r - 27r + 7r ( V T decays. The 
predictions of the Monte Carlo simulation are shown for the resonant and non-resonant r decays, and the background from 
other r decays and qq events. The resonant decays include decays with correct topology and a resonance (rj, fi or uj) in the 
final state. 



the same value [23[ . The result measured in this work is 
approximately 50% of the predicted rate but the ratio of 
the two branching fractions is consistent with unity. 



D. Non-resonant decay modes 

The resonant modes, involving 77, w and /1 mesons, 
do not account for all of the observed decays, as dis- 
cussed below. We consider the excess in the observed 



decays to be from "non-resonant" modes. We make no 
attempt to identify the contribution of resonances with 
larger widths (> 100MeV/c 2 ) as the nature of these res- 
onances is complex and their lineshapes will be modified 
by the limited phase space in the r decay. The Monte 
Carlo simulation describes the final states using a phase- 
space model for the final-state particles. The only excep- 
tion is the t~ — >• 37t - 27t + j/ t mode, which Tauola models 
using r~ — > a^v T decays [24j ]. 

We measure the branching fractions of the non- 
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TABLE V: Results and branching fractions for r — > 2tt iv + 3iv i> T , r — > Sir 2tv + u t , and r — > 3iv 2iy + tt°v t non-resonant 
decays. 



Decay Mode 


r" -> 2ir~ir + 3n°u T 


t — > 3tt 2n + i' T 


t" -> 3ir~2ir + n°u T 


Branching fractions (10 -4 ) 


0.10 ±0.08 ±0.30 


7.68 ± 0.04 ± 0.40 


0.36 ±0.03 ±0.09 


Data events 


4094 ± 64 


68985 ± 263 


7296 ± 85 


Efficiency 


(0.88 ±0.01)% 


(7.98 ± 0.02)% 


(3.71 ± 0.03)% 


Background 








Resonant 


1763 ± 222 


4441 ± 370 


4458 ± 244 


Other t decays 


1681 ± 44 


10621 ± 719 


1315 ± 100 


11 


573 ± 50 


1171 ± 205 


359 ± 37 


Total Background 


4017 ± 232 


16233 ± 835 


6132 ± 267 


Systematic uncertainties (%) 








Tracking efficiency 




3.8 


3.8 


7T° and rj PID 






3.0 


Pion PID 




2.5 


2.5 


Lepton-tag PID 




1.6 


1.6 


£ de- 




1.0 


1.0 


selection efficiency 




0.2 


0.9 


Background 




1.6 


22.9 


Total (%) 




5.2 


23.7 



resonant t~ — > 2tt^tt + 3tt°i' t , t~ — > 3ir~2ir + i' T , and 
t~ — > 3ir~ 2ir + ir° v T decays. The numbers of candidates 
are given by the numbers of events found in the data af- 
ter subtracting the resonant contributions and the back- 
ground from other r decays and qq events (see Table [Vj) . 

The invariant mass plots in Fig. [5] show that the res- 
onant decays dominate the r~ — > 2ir~it + ?)i{ a v T mode. 
The background is primarily from t~ — > ■n~it (jjv T and 
qq events. The branching fraction of the non-resonant 
t~ — > 27r _ 7r + 37r ( V T channel is determined to be (1.0 ± 
0.8 ± 3.0) x 10~ 5 . The systematic uncertainty on the 
branching fraction is dominated by the uncertainty in the 
background, which includes the Monte Carlo statistical 
uncertainty and the r branching fraction uncertainties. 
We do not list the fractional systematic uncertainties for 
this mode in Table fVl due to the smallness of the branch- 
ing fraction. The branching fraction is consistent with 
zero and we set a limit of 

B(t~ -> 27r-7r+37rV) < 5.8 x 10~ 5 

at the 90% confidence level. 

We also determine the inclusive r~ — > 2ir~ir + 3ir°i> T 
branching fraction, given by the sum of the resonant and 
non-resonant terms. We find the result (2.07 ± 0.18 ± 
0.37) x 10~ 4 , where the systematic uncertainty accounts 
for correlations between the systematic uncertainties of 
the individual modes. 

The t~ — > 2>-k~2-k + v t decay has only a small contribu- 
tion from resonant decays (see Fig. [B]). The branching 
fraction of the non-resonant r~ — > 2>-k~2-k + u t decay is 



determined to be 

B(t- ->• 37r-27r+^ r ) = (7.68 ± 0.04 ± 0.40) x 10~ 4 . 

The t~ — > 2Tt~ir + uji> T (a)—} 7r _ 7r + 7) decay is considered 
as a resonant background and is not included in the non- 
resonant branching fraction. Although the modeling of 
the 37r - 27r + invariant mass distribution is deficient, the 
selection efficiency remains the same if the Monte Carlo is 
re-weighted to resemble the data distribution. The decay 
model represents a significant improvement compared to 
a phase-space model, in which the p meson, observed in 
the 7r + 7r~ invariant mass spectrum, is not included. Fur- 
ther tuning of the model is required, which is outside 
the scope of the present study. The background from 
the qq events is validated by comparing the numbers of 
data and Monte Carlo events in the region above the 
t lepton mass. The inclusive t~ — > Z , n~2'K + v T branch- 
ing fraction is (8.33 ± 0.04 ± 0.43) x 10~ 4 and is obtained 
by adding the non-resonant branching fraction with the 
resonant branching fraction for the t~ — > tt~ f\v T via 
fi —> 2tt + 2tt~ decay. The branching fraction of the 
r~ — > ih~2h + v T decay (where h~ is either a 7r~ or K~) 
was measured to be (8.56 ± 0.05 ± 0.42) x 10~ 4 in a pre- 
vious BABAR analysis [l8| using a smaller data sample, 
which did not use charged particle identification and in 
which the branching fraction included the contribution 
of the t~ — > 2tt~ tt + ljv t decay. 

r~ — > 3ir~ 2ir + ir° v r decays are dominated by the res- 
onant modes (see Fig. [7]) . We determine the branching 
fraction of the non-resonant r~ — > 3n~ 2ir + ir a v T decay 
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FIG. 8: The n + ir~r] invariant mass in t~ — s> 27r~7r + 77^ T de- 
cay candidates for the 77 — > 77, 77 — > iv + it~iv , and 77 — ► Sir 
decay modes after all selection criteria are applied. The fit to 
the 77' peak (in the top two plots) is represented by the solid 
line. The number of 77' candidates in the 77— » Sir channel is 
determined by counting the number of events in a single bin 
centered on the 77' mass and subtracting the combinatorial 
events. The number of combinatorial events is determined by 
fitting the data (the solid line in the bottom plot), exclud- 
ing the data point near the rf peak. The peak in this plot 
indicates the expected location and width of an 7/ signal. 



mode to be 

B(t~ -> 37r-27r+7rV) = (3.6 ± 0.3 ± 0.9) x lO" 5 . 

The systematic uncertainty on the non-resonant — > 
3ir~ 2n + it v T branching fraction is dominated by the 
large uncertainty in the background (see Table |V| . Al- 
though the invariant mass distributions of the resonant 
modes in the Monte Carlo simulation are corrected to 
provide better agreement with the data, the corrections 
make little difference to the final branching fraction re- 
sult. The other r decays and the qq events contribute to 
a lesser extent; their contribution to the uncertainty of 
the background is very small. 

The t~ — > 37r~ 2ir + ir° 'v T (including 10 and excluding 
rf) branching fraction is (1.11 ± 0.04 ± 0.09) x 10~ 4 and 
is obtained by adding the non-resonant branching frac- 
tion and the resonant branching fraction attributed to 
the t~ — > 2ir~ir + u>h' T via u — > 7r~7r + 7r° decay. In ad- 
dition, the inclusive r~ branching frac- 



tion is (1.65 ± 0.05 ± 0.09) x 10" 4 and is similarly ob- 
tained by adding the non-resonant branching fraction 
and the resonant branching fractions attributed to the 
t~ — > 27T~7T + ?7J7 T via 77— > 7r + 7T~7r° and r~ — > 2ir~ir + Lui' T 
via uj— > 7r _ 7r + 7r° decays. 

The r~ — > 3ir~ 2ir + 7r° v t (including lu and excluding rf) 
branching fraction can be compared with isospin model 
predictions [III [26j]. There are three r decay modes 
with six pions in the final state: r~ — > 27r~7r + 37r c V r , 
t~ — > 3tt~ 2w + it v T , and r~ — > ■K~5if > v T (there are no 
measurements of the r~ — > n~5ir v T decay mode). There 
are four possible isospin states for six pion decays: (47rp), 
(3/?), (37T£j), and (irpuj). The relative rates of the de- 
cays can be used to identify the dominant isospin states. 
The approximate equality of the t~ — > 27r~7r + 37r°?y T and 
t~ — > 37r _ 27r + 7r°7y r branching fractions suggest that the 
(47rp) and (ftpuj) modes should dominate. The limited 
phase space imposed by the r mass suppresses the higher 
mass states and as a result we do not observe evidence 
of the p meson in these decays. 



E. Search for decays involving rf (958) decays 



We also search for the r 



7r-7rV(958K, 



-?- 

K~r]' (958) v T , and t~ — > 7r~?7'(958)^ r decays, where 
rf — > Tr~7r + r]. The first two decays are allowed first-class- 
current decays whereas the last decay is a second-class- 
current decay, with a rate that would be zero in the limit 
of perfect isospin symmetry. 

The numbers of rf candidates in the data and back- 
ground Monte Carlo samples are given in Table IVT1 For 
the t" — > 7r~7r°?7' (958)^ r via 77 — > 77 and the — > 
ir~ rf (958)zv r via r\ —¥ 77 and 77 — > 7r + 7r~7r° modes, we 
measure the number of 77' candidates by fitting the peak 
with a Gaussian function and the combinatoric back- 
ground with a polynomial function. The number of 77' 
candidates in the other channels is determined by count- 
ing the number of events in a single bin centered on the 
rj' mass and subtracting the combinatorial events. The 
level of the combinatorial background is estimated by fit- 
ting the mass spectrum or from the average level of the 
sideband region around the 77' peak. 

The 7r + 7r~?7 invariant mass distributions for the t~ — > 
-K~rf (958)zv T candidate decays are shown in Fig. [5] Al- 
though we see an rf peak in the r\ — > 7r + 7r~7r° channel, 
we find that it can be fully accounted for by qq events. 
We do not show the invariant mass distributions for the 
t~ -> 7r~7r°?7'(958)j/ r and r~ -» K~r)'(958)i> T decays. 
The analysis of these decay modes uses only the 77 — > 77 
and 77 — > 7r + 7T~7r° channels. The 77 — > 3ir° channel was 
not considered due to the limited size of the samples. 

The results for the three decay modes are given in Ta- 
ble IVII The selection efficiencies are determined with the 
signal Monte Carlo samples. The backgrounds from 77' 



mesons are attributed to e + e 



qq events and esti- 



mated using the Monte Carlo samples. The background 
estimations are validated by comparing the prediction of 
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FIG. 9: The K~2n~2n + , K + 3ir~n + , K~K + 2tt~tv + , K~2n'2n + n', and K + 3n~TT + n" invariant mass distributions in the 
data sample after all selection criteria are applied. The unshaded histogram represents r decays in which a charged pion is 
misidentified as a charged kaon, and the shaded histograms are primarily from qq events in which there is a charged kaon in 
the final state. The Monte Carlo simulation does not include signal decays. 



the Monte Carlo simulation with data for events where 
the invariant mass of all the observed final-state particles 
is greater than the r mass. 

We find no evidence for t~ — > 7r _ 7r°7/(958)i' T , r~ — > 
K~rf (958)^ T , or t~ — > n~r)' (958)v r decays (see Ta- 
ble IVip and place the following upper limits on the 
branching fractions at the 90% confidence level: 

B(t - -> 7r-7rV(958K) < 1.2 x 1CT 5 , 
B(r"-> K-r]'{958)is T ) < 2.4 x 1CT 6 , 
S(r - -)- 7r~?7'(958)^ r ) < 4.0 x 10~ 6 . 

The limits arc determined from the weighted average of 
the branching fractions measured for each mode. The 
t~ — > 7r~7r ?/(958)jA- and r~ — > K~i]' (958)is T channels 
are potential backgrounds to the r~ — » 7r~r/(958)^ r 
decay. We find that background from these two de- 
cays is less than two events based on the upper limits 
on the branching fractions and we consider these back- 
grounds to be negligible. The previous limits on the 
t~ — > w~ii'(958)u T decay were measured by BABAR to be 
7.2 x 10~ 6 Qj| and by CLEO to be 8 x 10" 5 [HJ. It is pre- 
dicted that the branching fraction of t~ — > tt~ r/' (958)i/ T 
should be less than 1.4 x 10" 6 [27| . 



F. Searches for decays involving charged kaons 

Finally we present the first search for high-multiplicity 
r decays with one or two charged kaons. We find no 
evidence for signal decays and place upper limits on 
the branching fractions of the t~ — > K~2ir~2ir + i' T , 

T~ — > K + ?>TT~'K + V T , T~ — > K~ K + 2k~'K + V t , T~ 

K~2-k~2'k + 'k v t , t~ -> K + 3n~Tr + TT v r , and r~ -> 
K~T}' (958)^ T decay modes (the t~ — > K~rj' (958)v T decay 
was presented in Section E). 

The events are divided into topologies in which the 
charged kaon has either the same or opposite charge as 
the parent r lepton. If there are two kaon candidates, 
they must have opposite charge. All other tracks are re- 
quired to be identified as charged pions. The selection cri- 
teria and systematic uncertainties are described earlier. 
The requirement on the invariant mass (M < 1.8 GeV/c 2 ) 
of the final state uses the kaon mass for tracks identi- 
fied as charged kaons. Figure [9] shows the mass spectra 
for the various channels. The predictions of the Monte 
Carlo simulation are divided into decays with or with- 
out a K~ (in this latter case, a ir~ is misidentified as a 
K~). The figures do not include any signal decays in the 
Monte Carlo samples. The background estimates, which 
give the dominant systematic uncertainty, are verified by 
comparing the numbers of events in the data and Monte 
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TABLE VI: Results and branching fractions for r — > ty 7r°f/ (958) v T , r — > K r/(958)^ T , and r — > n r]'(95S)u T decays. 



T — r TV TV Tj \\)oO)U T 


7] — ¥ 7T 


7) — y TV^TV TV^ 




Limit (yU/o b.L.J 


i a w 1 n — 5 
1.4 X 1U 


i.y X 1U 




Branching traction (1U J 


7.8 ± 4.1 ± 1.7 


(J.U ± 7.0 ± 9.o 




Data events 


O/i 4- in 

Z4 It 1U 


± D 




Background events 


5±7 


5±8 




Selection efficiency 


( 1 _i_ n no^Q^, 

IX. Do ztz U.UZ J /o 


I J..UU ztz U.Uo J /o 




r —f K rj [\)oo)}s T 


j] — > *yy 


7} — > TV TV TV 




Limit (90% C.L.) 


2.4 x 10 


4.2 x 10 




Branching fraction (10 J 


0.5 ± 1.3 ± 0.4 


1.6 ± 1.4 ± 1.2 




Data events 


a _i_ ^7 


15 ± 4 




Background events 


3±4 


11 ±3 




Selection efficiency 


(3.47 ±0.03)% 


(3.09 ± 0.04)% 




T _ -> 7T-77'(958)^ T 


77 — > 77 


7y— > 7r + 7T - 7T 


?7-> 3tt° 


Limit (90% C.L.) 


5.2 x 10" 6 


9.0 x 10" 6 


1.3 x 10" 5 


Branching fraction (10 -6 ) 


-2.8 ±3.5 ±1.9 


-0.4 ±3.9 ±4.3 


-1.8 ±8.1 ±3.3 


Data events 


40 ±22 


44 ±11 


12 ± 10 


Background events 


58 ± 12 


45 ± 12 


14 ±4 


Selection Efficiency 


(4.06 ± 0.35)% 


(3.25 ±0.15)% 


(0.96 ±0.10)% 



Carlo samples in the M > 1.8 GeV/c 2 region. 

The numbers of events selected in the data and Monte 
Carlo simulations are given in Table I VIII The back- 
grounds predicted by the Monte Carlo simulations are 
approximately equal to the numbers of events in the 
data sample. There is an excess of data events in the 
t~ — > K~2TT~2Tr + ir°i' T mode, but this excess extends to 
mass values above the r mass, indicating that events are 
due to background r decays or qq events. 

The upper limits on the branching fractions are given 
in Table IVIII There are no predictions for these modes. 
We estimate that B(t~ ->■ K-2n-2n + v r ) ~ KT 5 - 1CT 6 
if the decay is related to B(t~ — ? 3tt~2it + i' t ) by the ra- 
tio of the CKM matrix elements (V us /V u d)- The t~ — > 
3ir~ 2ir + ir a 'v T decay is dominated by decays to the nar- 
row low-lying resonances and the branching fraction of 
decay modes created by replacing a 7r~ with K~ would 
be highly suppressed due to the limited phase space. 



IV. SUMMARY 

We present measurements of the branching fractions 
for r lcpton decays to three-prong and five-prong final 
states. The results are summarized in Table IVIII1 The 
branching fractions exclude contributions of the K® me- 
son. The results are more precise than previous measure- 
ments and many decay modes are studied for the first 



time. 
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TABLE VII: Results and branching fractions for charged kaon decay modes. 



Decay Mode r 


"->■ K'2^2% + u T 


t — > K + 3tv ir + v T 


t~-> K-K + 2n-K + u T 


Limit (90% C.L.) 


2.4 x 10~ 6 


5.0 x 10~ 6 


4.5 x 10~ 7 


Branching fraction (10 -6 ) 


0.6 ± 0.5 ± 1.1 


1.6 ±0.6 ±2.4 


0.30 ±0.10 ±0.07 


Data events 


1328 ± 36 


1999 ± 45 


32 ±6 


Background 


1284 ± 72 


1890 ± 163 


15 ±4 


Selection Efficiency 


(7.9 ±0.1)% 


(7.9 ±0.1)% 


(6.7 ±0.1)% 


T~ 


-> K-2ir-2ir+w°u T 






Limit (90% C.L.) 


1.9 x 10~ 6 


8 x 10" 7 




Branching fraction (10 -6 ) 


1.1 ±0.4 ±0.4 


-0.6 ±0.5 ±0.7 




Data events 


112 ± 11 


154 ± 12 




Background 


84 ± 10 


170 ± 16 




Selection Efficiency 


(2.9 ±0.1)% 


(2.9 ±0.1)% 
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TABLE VIII: Summary of branching fractions excluding contributions from K° — > tv + tv . 



Decay Mode 


Branching fraction 


Resonant decays 






T — ^ 97T _ TV^ T)1J~- ritlplllHlTIP' /*1 ^ 

i 7 — m J i 1 lis f I lllLil l_l 111 t i- I 


(2.25 ±0.07 ±0.12) x 10" 


-4 


t~ — )■ i l'K 'R Jr T\v T (excluding /i) 


(0.99 ±0.09 ±0.13) x 10" 


-4 


t~ — ^ tt ^Tr^/ykv (including /i) 


(2.01 ±0.34 ±0.22) x 10" 


-4 


r~ — > tt~ fiv-r via /i — > 27r + 27r~ 


(5.20 ±0.31 ±0.37) x 10" 


-5 


r~ — > ir~ fii/ T via /i — > TV + ir~r) 


(1.26 ±0.06 ±0.06) x 10" 


-4 


B{fl^27V + 27V-)/B{fl -> TTTVr/) 


0.28 ±0.02 ±0.02 


T~ — i> 27T~ TV^LdVr 


(8.4 ± 0.4 ± 0.6) x 10" 


-5 


t~ -> 7r~27r°wz/ T 


(7.3 ± 1.2 ± 1.2) x 10" 


-5 


Non-resonant decays 






t~ — ^ 37r ~27r^~z>' T (excluding cj, ^i) 


(7.68 ±0.04 ±0.40) x 10" 


-4 


r~ — > 27r _ 7r + 37r°i/T- (excluding 77, cj, /1) 


(1.0 ±0.8 ±3.0) x 10" 


-5 


r~— > 27r _ 7r + 37r°!/ T (excluding 77, /1) 


(16.9 ±0.8 ±4.3) x 10" 


-5 


r~— > 37r _ 27r + 7r t V T (excluding 77, cj, /1) 


(3.6 ±0.3 ±0.9) x 10" 


-5 


r~— > 37r _ 27T + 7r°!/ T (excluding 77, /1) 


(1.11 ±0.04 ±0.09) x 10" 


-4 


Inclusive decays (including 77, a;, /1) 






r~ — > 2n~iv + 3n v T 


(2.07 ±0.18 ±0.37) x 10" 


-4 


t~ — > 37r _ 27r + i' T (excluding oj) 


(8.33 ±0.04 ±0.43) x 10" 


-4 


r~ 3tt~ 2tt + it v T 


(1.65 ±0.05 ±0.09) x 10" 


-4 


t! (Q^r\$\ HppfivQ iiT^T^pr IpvpI pnnfidpnpp limits 

11 i l/iju 1 ucLa y id y-yKJ/v lij^jjcl icvci l.uiiiiuciil c 111111 ii i 






t~ -nr-n°r)'{958)u T 


< 1.2 x 10" 


-5 


r"-> Jf"r?'(958)^ T 


< 2.4 x 10" 


-6 


T--S- 7T-?7'(958)^ T 


< 4.0 x 10" 


-6 


Kaonic decays (90% upper level confidence limit) 






t~ — ^ K-2n~2ir + v T 


< 2.4 x 10" 


-6 


t~ — > K + 3-k~tt + v t 


< 5.0 x 10" 


-6 


T~ -> iY _ Jsr + 27T _ 7r + !/ T 


< 4.5 x 10" 


-7 




< 1.9 x 10" 


-6 


r _ -> K+37r-7r+7r ( V T 


< 8 x 10" 


-7 



